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Abstract 



Nonthermal fixed points represent basic properties of quantum field theories, 
in addition to vacuum or thermal equilibrium fixed points. The functional 
renormalization group on a closed real-time path provides a common 
framework for their description. For the example of an 0{N) symmetric 
scalar theory it reveals a hierarchy of fixed point solutions, with increasing 
complexity from vacuum and thermal equilibrium to nonequilibrium. 



1 Introduction 



Thermal equilibrium properties of many-body or field theories are known 
to be efficiently classified in terms of renormalization group fixed points. A 
particularly powerful concept is the notion of infrared fixed points, which 
are characterized by universality. These correspond to critical phenomena 
in thermal equilibrium, where a characteristic large correlation length leads 
to independence of long-distance properties from details of the underlying 
microscopic theory. In contrast, a classification of properties of theories 
far from thermal equilibrium in terms of renormalization group fixed points 
is much less developed. The notion of universality or critical phenomena 
far from equilibrium is to a large extent unexplored, in particular, in 
relativistic quantum field theories. Here the strong interest is mainly driven 
by theoretical and experimental advances in our understanding of early 
universe cosmology as well as relativistic collision experiments of heavy nuclei 
in the laboratory. 

It has been recently demonstrated that nonthermal fixed points 
exist in relativistic self-interacting scalar field theories, and that the 
associated far-from-equilibrium long-distance properties are characterized 
by universality [1_J. Compared to a thermal distribution at low momenta 
~ p~^, a nonthermal scaling solution was shown to exhibit strongly enhanced 
fluctuations ~ in the infrared for three spatial dimensions. From a 
large class of initial conditions leading to nonequilibrium instabilities, such 
as spinodal decomposition or parametric resonance, these nonthermal fixed 
points are approached without fine-tuning of parameters. As a consequence, 
they may have an important impact on properties of early universe dynamics 
after inflation, or on the physics of thermalization in heavy-ion collisions. 

In this work we employ the functional renormalization group for the 
average action [2] to determine properties of nonthermal fixed points 
for 0{N) symmetric scalar quantum field theories. This provides a 
complementary approach to the two-particle irreducible (2PI) effective action 
techniques that were used in Ref. The functional renormalization group 
has a long record of successful applications to the nonperturbative physics of 
critical phenomena in thermal equilibrium as reviewed in Refs. [31 H]. While 
in equilibrium typically a Euclidean formulation is adequate, nonequilibrium 
properties require real-time descriptions. For quantum systems a generating 
functional for nonequilibrium correlation functions with given density matrix 
at initial time can be written down using the Schwinger-Keldysh closed 
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time path contour following standard procedures [5]. In principle, this 
can be used to construct a nonequilibrium functional renormalization group 
along similar lines as for Euclidean field theories in thermal equilibrium. 
However, important differences include the absence of time-translation 
invariance for general out-of-equilibrium situations [6]. The nonequilibrium 
renormalization group takes on a particularly simple form at a nonthermal 
fixed point, since it becomes independent of the initial density matrix and 
we concentrate on this in the following. 

Far-from-equilibrium critical phenomena can be observed in many 
systems, ranging from self-organized criticality in the well-known Bak-Tang- 
Wiesenfeld sandpile model [7| to phase transitions between nonequilibrium 
stationary states for species dynamics in biologically motivated descriptions 
|8]. Renormalization group techniques have been employed to nonequilibrium 
critical dynamics in a variety of formulations, with a large body of 
work concerning applications to classical react ion- diffusion problems [9l 
[To] , stationary-state transport through quantum systems coupled to 
reservoirs [HI [T21 [131 [E], or cosmology PSJ. 

Our paper is organized as follows. Sec. 11.11 gives an introduction to 
nonthermal fixed points. Sec. [2] discusses the exact renormalization group 
on a closed time path, where 12.11 describes the connection to the 2PI 
effective action for later comparison. Though the derivation is in principle a 
straightforward generalization of well-known vacuum or thermal equilibrium 
techniques as pointed out in Refs. [IHl [131 HI], there are some important 
differences. This concerns the exact renormalization group equation for 
the scale-dependent generating functional and its relation to standard 
diagrammatic expressions. From the functional form of the renormalization 
group equation given in Sec. 12.21 the flow of all correlation functions can 
be obtained by further functional differentiation. In contrast, the standard 
nonequilibrium diagrammatic rules derived in Sec. l2.3l exploit the vanishing of 
anomalous propagators and vertices at extrema of the generating functional. 
Another important difference is the apparently larger freedom of choosing 
cutoff functions on a closed time path. Sec. 12.21 discusses constraints which 
arise from constructing a proper renormalization group flow that interpolates 
between short and long distance scales. The analysis of nonthermal fixed 
points is given in Sec. [31 Sec. 13.11 describes the scaling form of correlation 
functions, which are evaluated from the flow equations using a resummed 
large-iV expansion to next-to-leading order (NLO) in Sec. 13.21 The fixed 
point results are discussed in Sec. 13.31 Conclusions are given in Sec. [H Two 
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Appendices discuss details about the functional form of propagators and a 
specific representation for cutoff functions. 

1.1 Thermal versus nonthermal fixed points 

Important questions about the evolution of the post-inflationary early 
universe or heavy-ion collision experiments concern thermalization processes. 
Starting from a given state or density matrix at some initial time, one 
considers the nonequilibrium time evolution towards thermal equilibrium. 
The time to approach the latter can crucially depend on the presence of 
infrared fixed points at which time scales diverge. This is a well known 
phenomenon near thermal second-order phase transitions where the dynamics 
is characterized by critical slowing down. Nonthermal fixed points can 
have the dramatic consequence that a diverging time scale exists far from 
equilibrium which prevents or substantially delays thermalization, depending 
on how closely the fixed point is approached during the nonequilibrium time 
evolution. 

Fixed points correspond to time and space translation invariant scaling 
solutions for correlation functions. Here we consider a relativistic scalar 
field theory with Heisenberg field operator For instance, the two-point 

correlation function G{x,y) can be identically written as 

G{x,y) = ($(x)<l>(y))e(x° + ($(y)$(x))e(y° 

= F {x, y) - -p{x,y) sign{x^ -y^) , (1) 

where (...) involves the trace over a given density matrix and x = {x^,x.) 
with time and vector x in c? spatial dimensions [5]. The statistical two- 
point function F{x, y) determines the anti-commutator, while the spectral 
function p{x, y) describes the commutator of two fields according to 

F(x,i/) = ^({$(x),$(y)}) , p{x,y)=imx),^y)]). (2) 

Loosely speaking, the spectral function p(x, y) determines which states are 
available, while the statistical propagator F{x, y) contains the information 
about how often a state is occupied. 

A tremendous simplification of thermal equilibrium is that the spectral 
and statistical two-point functions are related by the fluctuation-dissipation 
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relation. In Fourier space with momentum p the latter reads 

-^(P) I thermal equilibrium = ^^t(P°) + p{p)\ thermal equilibrium ; (3) 

where nrip^) = (exp(p°/T) — 1)^^ denotes the Bose-Einstein distribution for 
a temperature T. In the infrared this distribution is given by 

In this work, we consider scaling behavior at an infrared fixed point that 
can be expressed in terms of the standard dynamical scaling exponent z, the 
anomalous dimension t] and an "occupation number exponent" k, according 
tcfl 

F(/,p) = .2+-F(sy,sp) , p{p',p) = s'-'^p{sy,sp). (5) 
Inserting ([5]) into ([3]) leads to the following relations for 

vacuum (T = 0) : k = —r] , (6) 
thermal equilibrium (T ^ 0) : k, = —i] + z . (7) 

In contrast, out of equilibrium F{p) and p{p) are not constrained by a 
fluctuation-dissipation relation. This additional freedom allows new scaling 
solutions, which do not exist in thermal equilibrium because of the presence 
of ([3]). We show in Sec. [3] that one can find at a 

nonthermal fixed point : k, = —rj + z + d , (8) 

using nonperturbative large- iV expansion techniques to NLO. In contrast to 
the equilibrium relations ([6]) and ([7]), out of equilibrium the dimensionality 
of space d enters in ([8]). This leads to comparably large values of 
which corresponds to strongly enhanced infrared fluctuations as observed 
in Ref. [T]. As for the thermal equilibrium case, we show that nonthermal 
scaling solutions correspond to fixed points of the functional renormalization 
group. In order to observe both thermal and nonthermal fixed points the 
functional renormalization group has to be considered on a closed time path, 
which implements the possible absence of a fluctuation-dissipation relation. 



-'^In Ref. [T] this exponent is called a following literature on ultraviolet (UV) fixed points 
related to turbulence (see also Sec. l3.3[) . Since here we concentrate on infrared fixed points 
associated to critical phenomena we use k, which cannot be confused with the specific heat 
exponent. 
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2 Functional RG on a closed time path 



2.1 Generating functional 

We consider a scalar A^-component quantum field theory with 0{N) 
symmetry, whose classical action reads 

S[p\ = \j ^^a{x)iD-^{x,y)^b{y) - ^ j ^^a{x)Va{x)^b{x)Vb{x). (9) 

Here the free classical inverse propagator is given by 

iD-2{x, y) = - + m") 5ab5^'''-'\x - y) , (10) 

with Dx = d^o—d^- Summation over repeated indices a = 1, . . . , is implied 
and we use the shorthand notation ^ = da;" J d'^x. For nonequilibrium 
initial-value problems, the closed time path C leads from the initial time to 
along the real time axis and back to to [16j. For the discussion of stationary, 
i.e. time translation invariant nonequilibrium solutions to ~^ — oo, which we 
will consider in the following. 

We follow closely Refs. [ITJ [18] and employ a standard notation, where 
the superscripts '+' and '— ' indicate that the fields are taken on the forward 
(C+) and backward (C^) branch of the closed time path, respectively. Then 
the action can be written as 



A 
4!iV 



{'^aix)vUx)vtix)'^tix) - Va i.x)Va i.^)Vb i.^)Vb i.^)) . 

(11) 



where the minus sign in front of the '— ' terms accounts for the reversed 
time integration with = dx° J d'^x = J d^^^x. The corresponding 
generating functional for correlation functions then reads: 

Z[J,R] = I Dv^+P^-expjz S[v\v-]+l^ (^«(^)'^-~^^)) (-}-%) 

^b (y) 



(12) 
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Here we introduced apart from the source terms hnear in the fields, 
and J~, also sources bihnear in the fields, etc. Again the superscripts 
'+' and ' — ' indicate that the respective time arguments of the sources are 
taken on the forward or backward branch of the closed time path. As a 
consequence of this notation all space-time integrations run from — oo to oo, 
thus calculations are done in d + 1-dimensional Minkowski space-time. 

The connection to statistical and spectral components as employed in 
Sec. 11.11 is conveniently obtained with a linear transformation A of the fields 

-a('C). (13) 



where 



We note that (p agrees with the defining field in iQ only for ip~^ = (p~ . Since 
this will be the case for expectation values in the absence of sources and 
since there is no danger of confusion in the following we keep this notation. 
Correspondingly, we write for the source terms 

rv^. (15) 



J J \ J- J ' \ R'^ J \R-^ R 
Inserting these definition^ into the functional integral f[T^ yields 



Ja{x) 
Ja{x) 



+ 



(16) 



^Note that (fT5|) can be equivalently written as 
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where S[ip,ip\ = So[ip,ip] + Si^t[f,f] consists of the action for the free field 
theory 

(17) 

and the interaction part 

Smt[^,^] = J <paix)iPaix)ipbix)ipb{x) - J ipaix)ipaix)ipbix)ipbix) . 

(18) 

A similar construction can be done for the corresponding classical-statistical 
field theory, i.e. neglecting quantum corrections. For a detailed derivation 
see Refs. [T71 IT9] . The result for the classical-statistical generating functional 
is given by (|T6l) and f|T7|) with S'int replaced by 



Sint[^,<P] = j^^a{x)^a{x)^b{x)^b{x) • (19) 

Comparing 0181) and f[T^ one observes that the quantum theory includes the 
classical interaction vertex, however, in addition there is a vertex present 
which is absent classically. 

From the generating functional for connected correlation functions 

W = -i\nZ (20) 

we define the macroscopic field 0a, and (pa by 

-J— = 0a(x) , — — = 0,(a;). 21 

The connected retarded/advanced propagators G^^^{x,y)^ the statistical 
correlation function Fab{x,y), and Fab{x,y) are defined by 



dJa{x)6Jb{y) 6Ja{x)6Jb{y) 

,j , , . = iFabix,y) , = zFabix,y), (22) 

6Ja{x)6Jb{y) 6Ja{x)6Jb{y) 

and one observes 

Gtbix,y) = Gliy,x). (23) 
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These are related to the spectral function defined in ([2]) by 



pab{x, y) = C^t^ix, y) - G^t,{x, y) . 



(24) 



Similarly, we employ 



R. 



(25) 



We emphasize that the macroscopic fields and correlation functions 
depend on the sources. In particular, for vanishing sources J, J, 
R^, R^ the field = 0, the function F = 0, and, more generally, 
S'^W/SJ{xi)SJ{x2) ■ ■ ■ SJ{xn)\j^j^jiF=jiF^Q = llHj. It is important to note 
that this holds in the presence of non- vanishing retarded/advanced sources 
i?^'^, which can be observed from the fact that these can be subsumed 
according to f[T^ and flTTl) into a modified inverse classical propagator as 
iD~^ + R^'^. This will be used for the evaluation of the functional fiow 
equation in Sec. 12.21 

The Legendre transform of W[J, J,R^'^'^'^] with respect to the source 
terms defines the two-particle irreducible (2PI) effective action [201 [HI l2T] 



The 2PI effective action has been extensively used to describe nonequilibrium 
dynamics. For a review and generalization to nPI functional integral 
techniques in this context see Ref. [5]. A powerful nonperturbative 
approximation is the 2PI expansion to NLO [22l [23] . which 

was successfully applied to far-from-equilibrium dynamics and late-time 
thermalization [22l [211 ISSl [261 [I] as well as critical phenomena in thermal 
equilibrium |27]. In Sec. 13.21 we will recover the 2PI expansion from the 
functional renormalization group (see also Ref. [6J). The close connection 
between nPI and functional renormalization group techniques [4] will be 
exploited in the following for the discussion of nonthermal fixed points. 



^2PI[0,0,G^G 



= W - I Ua{x)Ja{x)+^a{x)Ja{x) 



J X 




(26) 
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2.2 Exact flow equation 

We consider the average action [2], evaluated along the closed time path as 
described above. The average action can be conveniently obtained from the 
2PI effective action fl26|) by partially undoing the Legendre transform with 
respect to the source terms quadratic in the fields: 

r [0, 0; ^Y2Pi-'-Tr[ G^i?^ + G^/?^ + tFR^ + zFi?^ } . (27) 

Here we employ a compact notation, where the trace involves Minkowski 
space-time or momentum integration as well as summation over field indices. 
The remaining dependence on the bilinear source terms R^'^'^^^ will be 
used in the following to construct the renormalization group flow of the one- 
particle irreducible (IPI) effective action. The equations of motion for the 
fields are 



-Jaix) - / (R^h{x,y)(j)biy) 

+ \Rab{^.y)Uy) + \Hy)RUy.< 



-Jaix) - / iR^h{x,y)(f)biy) 

+ l-Rlix, y)Uy) + lMy)Rtiy, ^)) , (28) 



6(f)a{x) 



and 

^Kbi^^y) 2 5Rab{.x,y) 2 

— = = -Fab(x,y) , — — — ^ — - = -Fnh(x,y) . 29 

5i?f,(x,i/) 2 '-'^ '^^ 5R^ab{^,y) 2 '^"^ '^^ ^ ' 

As a consequence of the partial Legendre transform the quantity defined 
by (|27|) becomes i?R''^'^''^-source dependent. For the renormalization group 
flow these source terms are taken to depend on some scale squared, k"^ > 0, 
which we denote as R^'^'^'^. The average action 

Tk[<PA] = r[0,0;i?i^'^'^'^] (30) 
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as well as all correlation functions are taken to be scale dependent functionals 
of the fields 4> and (p, such that F = 0] etc. The behavior under changes 
of the scale is described by the fc^-derivative of fl271) . With the notation for 
the dimensionless "beta function" derivative = 2k'^dTk/dk'^ this yields 

^k[<PA = -^Tr|G^[0,0]f?^ + G^[(/.,0]i?^ + zF,[0,0]i?f + |. 

(31) 

This exact flow equation has a standard one-loop form, where the scale- 
dependent loop-propagators are related to second derivatives of the average 
action, e.g. 

d(pa{x)d(pb{y) d(pa{x)d(pb{y) 



In the compact matrix notation employed above these relations are given by 
the identities: 





[(rf 


' + Rf) 


- ixf + Rk)ixt + Rt) 






\-' 


Gtl^A] = -\ 


[(rf 


' + Rt) 


-(rr + ^f)(rf + ^^) 




' + r!)] 


\-' 




[(rf 


' + r!) 


- {Vf + Rt){Tf + Rl) 




' + Rf)] 


\-] 




[(rf 


' + Rk) 


- iXT + Rt)iXf + Rt) 




' + Rt)] 





where we have used (1251) . These are obtained by noting that the relation 
between the fields and sources is fc-dependent, i.e. J = Jk[4>,(f>] and J = 
Jk[4>, 0], and a derivation is given in Appendix lAl 

The sources i?^'^'^'^ are typically chosen from a family of functions to 
construct a renormalization group flow that interpolates between the classical 
action S[(t>, cp] for fc^ cxd, or some UV cutoff scale fc^ ^ A^, and the standard 
IPI effective action r[0, 0] = r[0, 0; ^ g] for P 0: 

hm rfc[0,0] = 5[0,0] , limrfc[0,0]=r[0,0]. (34) 

This is standard for Euclidean field theories, however, on a closed time path 
there are few important differences which help to constrain the freedom in 
choosing R^''^'^'^ , which is discussed in the following. 
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Using the equations of motion for the fields (125]) with 



and the defining functional integral ( fT6ll with ( l20| ) and ( !26ll one can rewrite 
(EID functional integro-differential equation! 



J3 



S[(j), (f)] — i\n J VipVip exp yi 



S[(t) + ^,(j) + Lp] - S[(t), 



With the representation of the delta functional 



.(36) 



j V(f>expi^i j Lpa{x)(pa{x)^ = 5[(f] (37) 

one observes from (1361) that the interpolation property can be efficiently 
implemented with the class of functions 

R^abi^^ y) = Rk{-^x)5ab5{x - y) , i?fffe(x, y) = 0, (38) 

where in Fourier space for fixed momentum the infrared cutoff function is 
required to behave as 

lim i?fc(p^) ^ oo , lim i?fe(p2) ^ 0. (39) 

We emphasize that i?^' cannot be chosen to diverge as i?^' for large k 
since the latter has to dominate for short-distance scales in order to recover 
the 5-function behavior fl57|) . For Rj^'^ = one obtains from flHT]) 

r,[0,0] = -^Tr{G^[0,0]f2^ + Gt[0,0]i?^}. (40) 

■^For the derivation one uses that the measures are invariant under the shifts ip —> (j> + ip 
and ip —>■ (j) + (p. 
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This exact flow equation will be used below for the description of nonthermal 
fixed points. With 0231) and (1251) it can be written in a more compact 

form using Tr jc^j^, | = Tr |g'^[0, 0]_R^|. However, we will keep 

the explicit symmetry between retarded and advanced functions in the 
notation to simplify the diagrammatics below. The fact that in (140|) the 
"off-diagonal" retarded/advanced functions R^'^ are sufficient to implement 
an infrared cutoff on fluctuations coincides with the property that in the 
defining functional integral (ITB]) the free propagator (1171) exhibits vanishing 
"diagonal" entries. This motivates similar constructions of fiow equations by 
modifying the frequency or momentum dependence of free propagators as in 

Refs. nani]. 

Exact fiow equations for ra-point functions are obtained from n functional 
derivatives of fHUl) with respect to the fields and 0. For instance, 

ffi(-.y) = -^4^^^^fl? + ^^^^fl.4, (41) 

2 \^5(j)aix)5(f)biy) 6(j)aix)6(j)biy) J 

where the derivatives will be evaluated in Sec. 12.31 using (|33l) . Since the 
latter equations for G^'^lcp, (p] contain already second derivatives of Tk[4>,4>], 
one recovers the fact that fiow equations for n-point functions depend on up 
to {n + 2)-point functions. 



2.3 Diagrammatics 

After taking derivatives of (HOl) to obtair0 exact fiow equations for ra-point 
functions, all sources can be set to zero except for R^'^, which implement 
the scale dependence as discussed above. Following Sec. I2.1l this corresponds 
to = and = 0, while the latter entails Vf,'^ = according to 
(133]). More generally, 5'T[0, 0]/50(xi)50(x2) ■ ■ ■ 50(a;„)|^^^^ = 0, where 
the macroscopic field in the absence of sources, (p = (po, may be nonzero 
corresponding to spontaneous symmetry breaking [18]. From ( !33|) one finds 

GfJ4> = ^o,4> = o] = -(rt* + i?«)"' [0 = 00, = 0] = ^..A, 

\ / ab 

GtM = 00, = 0] = - (rt^ + Rt)l^ [0 = 00, = 0] = ^ , 

^i^M40 = 0o,0 = O] = (G^rf G^)^j0 = 0o,0 = O] =^ ^.(42) 

^The diagrammatic results derived from ((40|) or from ([3T|) agree, of course. 
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For a diagrammatic description a solid line is associated to the statistical 
correlation function iFk, and the retarded (advanced) propagator (G^) 
is represented by a solid-dashed (dashed-solid) line. The cutoff function 
derivatives and are associated to a cross, with reversed ordering of 
dashed to solid lines as compared to the retarded and advanced propagators: 

Rlab = ^--^ , RU = • (43) 

Similarly, the proper vertices obtained from n functional 0-derivatives and 
m functional ^-derivatives of Tk[4>, (p] are represented by full dots connected 
to n dashed lines and m solid lines, once the fields are evaluated at (f) = (po 
and (j) = 0. Four-point vertices are displayed, for instance, as 

c d c d 

k.abcd ' k.abcd • / ^ ^\ 

a-^b a- 'b (44) 

We will frequently denote proper four-point vertices as F^"'^^^^, irrespective of 

(2) 

their external lines, and similarly for two-point functions F^ etc. 

With this notation the exact renormalization group equation fHOj) reads 
for Ffc = Ffc[(/) = = 0] diagrammatically: 



(45) 



It should be emphasized that diagrammatic descriptions of flow equations 
employ ( H2l) . which exploits the vanishing of the anomalous propagator 
Fk and F^*^ at extrema of the generating functional. In contrast, the 
renormalization group equation fHOl) using the full functional form of the 
propagators fl33|) is still required to obtain the flow of ra-point correlation 
functions by further functional differentiation. 

In view of the application to nonthermal fixed points, we will consider 
in the following no spontaneous symmetry breaking (0o = 0) for simplicity, 
though the generalization to 0o 7^ is straightforward. In the absence of 
spontaneous symmetry breaking, where three-point vertices vanish for the 
theory described by (Q, we derive the diagrammatic expressions for F^ 

and F^^abcd froi^ derivatives of fBOj) employing the functional form of the 
propagators fl33|) . This is used to observe standard diagrammatic rules, which 
facilitate the computation of flow equations for higher n-point correlation 
functions. 
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According to ( 14TI) the exact flow equation for tlie two-point function 
involves the functional derivative 

^ ^k,ab[4>j 4>] _ J^R -p4><f>4>4> _ -pH4>4> fr't'i' J- -p4>4> /^R 

^k,ae'- k,efcd^kjb ^ k,ae'- k,efcd ~^ ^k,gh^k,hb 



-1 , 7 7 7 / - . N -1 



R 

k,ij'^k,jb 



- Gl^jtif (r*^^ + i?^) Tt%Gl,, } '^l ' (46) 

suppressing space-time arguments and neglecting three-point vertices. This 
is obtained from (l33l) with R^'^ = 0. Since Tff'[(j) = 0, = 0] =0, only the 
first and last term of (H6l) contribute to the compact diagrammatic expression 

for r^'^fe. Accordingly, for f = f [0 = 0,0 = 0] we find 

a b ah ah ah \^') 

Since a four-point function involves two further functional derivatives, it 
receives also contributions from derivatives of the anomalous terms r^'^[0, 0]. 
Evaluated at = = these do not vanish in general, since only 
Y^H^^ = 0, = 0] = 0. We obtain 

Z'^*** /■'X** z^^"* z**^'* z"^'* 



a h c a a h c a a h c a a h c a a h c a a h c a 

+ 



2 



i'X*^ ^'X^ *'^^*\ 

h"f\ ^ h-f\ ^ hrf\ ^ h"T\ ^ ^■r\ ^ /^rT\ 

a h c d a h c d a h c d a h c d a h c d a h c d 
+ P{a, h, c, d) I 

.>s^ 



/TVs, , , 

ah c d ah c d ah c d ah c d (48) 



where P(a, 6, c, d) denotes permutations of the respective field indices. 
Likewise, flow equations for six-point or higher vertices can be derived by 
further functional differentiation of (l40l). 
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Specific two- or four-point functions are then obtained from (1471) and 
by attaching the respective external dashed (0) or sohd (0) fines. For 
instance, starting from f[^Tl) the corresponding fiow equation reads according 
to fH7j) in Fourier space 

HUp) = ~U \ rKl(p,-p,?)G't/(?)^?/c(g)Gt.(g) 



+ rgtl(P, -P, <lYFu,ef{q)RifM)GtM \- (49) 



wfiere = j &^+^q/ {2TtY+^ . 

Tfie fiow equations can be simpfified considerably using specific 
representations of the cutoff functions R^'^ ■, such as a sharp cutoff discussed 
in Appendix [Bl There we show how the latter can be employed to integrate 
the flow analytically. This will be, in particular, important for more 
sophisticated approximations than the one we employ in Sec. [3l and we give 
the respective formulae in the Appendix for completeness. For our case, it 
is advantageous to proceed without restricting the calculations to a specifc 
cutoff and to note that (HUI) can also be written, using ([21]), as 



= Tr|flf^ + fl*^jrt|,#,,,^l = 4r»[,#,,,^]. (50) 

For instance, the r.h.s. of the flow equation for the two-point function (1471) 
can be rewritten using relations such as 

{Vf + Rf)~' R^ (rf + Rf)~' = -d, (Vt" + Rf) . (51) 



^We employ 

/ \ p n 



i=l 
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Similar relations are typically employed in a reversed sense, i.e. to derive 
renormalization group equations for n-point functions using diagrammatic 
one-loop techniques [3], SI |10]. In order to simplify the notation, we use 
GI'^ = Gf^[(l) = 0,0 = 0] and iFk = = 0,4> = 0], such that the 

retarded/advanced and statistical propagators denote the expressions P2|) . 
For instance, the equation for advanced functions corresponding to fIFI]) reads 
then 

G^R^G^ = dkG^ , (52) 

and we have 

GfRftFk + tFkRtGt = 4 ^Fk • (53) 
The flow equation for the four-point function fHSl) involves 

^k^k^k'- k ^k'- k + ^k'- k ^k^k^k'- k -(^ky^k'-k^k'-k ) ■> y^^) 
and 

^RdR^Rp(4) p(4) R (4)^RpR .p p(4) „R (4) dA^A 

'^k^k'^k^k ^-^ki-k +'^k^k '^k^k^^ki-k +<-^fcJ-fc ^^krik'^k 

= ^k{GlTf^FkTf), (55) 



as well as 



(GlRltFk + ^F,rf + zF.rf ^F, + iFkR^Gi 

= dk fzF.rf zF^rl")) . (56) 



Applying (15Ti) - (l5^ the flow equation for the two-point function (H71) is 
given by 

e = 44 {Q + Q + Q}. 

a 6 a 6 a b (57) 

Like before, equations for specific two-point functions can be directly deduced 
from ( 1571) by attaching the respective external dashed (0) or solid (0) lines. 
For instance, in Fourier space the two-point flow equation (ll9l) now reads 

^iUp) = -\dkj { Tt%{p,-p,q)Gl,M) 

+ 4%iv.-vA)GUM) 

+ rti'.(p,-p,?)^^Mc(g)j. (58) 
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The four-point vertex flow equation (H5]l with yields 



b'^'^d b'^^-^d 



P{a, b, c, d) 




2 

Trf (59) 

Again, we display all combinations of retarded, advanced and statistical 
propagators separately for illustration, though several obvious graphs are 
identical or vanish and ( l59l) can be given a more compact form. This is 
discussed in more detail in Sec. [3], where an approximation in terms of a 
resummed 1/A^ expansion to NLO is employed. 

The above equations for F^^^ and t^^^ have to be integrated from the 
ultraviolet scale A, with Tk=A = S, to any given infrared scale k, and r^^o = 
r as discussed in Sec. 12. 2[ For the quantum theory with action described 
by ( |T7I) and ( |T8l) the only nonvanishing two- and four-point functions at the 
ultraviolet scale are 

rt,(a:, y) = Tf;jx, y) = tD-l{x, y), (60) 

with the classical inverse propagator (|TOl) . and interaction vertices 

(x-w), (61) 

rSt.(^,2/,^,^) = \^i%{^,y,^,y^)- (62) 

Following the discussion of Sec. 12.11 we note that for the corresponding 
classical-statistical field theory the "quantum vertex" Y'^'^'^ vanishes 
according to (fTOl) . In contrast, in the quantum theory it is only suppressed 

by a factor of four compared to F^"^*^"^. The latter is present in both quantum- 
statistical and classical-statistical field theories [T7|n[9]. 
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3 Nonthermal fixed points 



3.1 Stationarity condition 

Infrared fixed points correspond to low-momentum scaling solutions for n- 
point correlation functions. As mentioned in Sec. II. ![ the absence of a 
fluctuation-dissipation relation out of equilibrium can lead to new nonthermal 
scaling solutions. This was pointed out recently in Ref. pQ for the theory 
considered here, based on a stationarity condition for time translation 
invariant solutions of nonequilibrium time evolution equations [5]. An 
equivalent (scale-dependent) stationarity condition can be obtained from the 
functional renormalization group, where it appears as a nontrivial identity 
relating the various second functional derivatives of rfc[0,0]. For this we 
consider 

^^fc^c {^fc^ca Gk,ca} = '~'^^f,ac^t'!'cd^k,de | (^fc ) ea ~ )ea } 

= i^Mc{rg,-rg,} , (63) 

where we suppress space-time dependencies in the notation. The first line 
of fl63p can be verified by simple matrix multiplication, and the second line 

follows from fl55]) and fH2]) . Following conventional notations, iVff' represents 

the statistical part and F^"^ — Vf^ the spectral part of the self-energy denoted 

as [m 

yiF _ ,n</'<^ yip _ T^<f><f> _ p0<^ (RA\ 

^k,ab — ''^ k,ab ' '^k,ab — ^ k,ab ^ k,ab ■ \^^) 

The latter is related to the retarded/advanced self-energy by = S^^^j, — 
S^^^, similar to the spectral function being related to the retarded/advanced 
propagator by pk^ab = G^ab ~ ^tab according to (^^. Using this notation 
( |63|) reads in Fourier space 

This equation is well-known in nonequilibrium physics. In the language of 
Boltzmann dynamics it states that "gain terms" equal "loss terms" for which 
stationarity is achieved [S]. 

The flow equations for the self-energies are determined with flM|) and fl57|l : 

Kabip) = l^kj^ {rgt(p, -P, q)GlM + rfSi(p, -P, q)Gt,M 
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+rrr.(p,-p,?)^^Mc(g)|, (66) 

rS.(p, -P, q) - rttt(p, -p, q)] Gl,M) 



+ 
+ 



fc,a6cd(P' -V. q) - ^labM -p^ q) G.dM 



KabM -p^ q) - ^k abM -p^ q) iFk,dciq) 



(67) 



For given four-point functions the stationarity equation ( l65l) with ( l66l) and 
del!) depends on Fk,ab{p) and G^„fc(p) = G^^,^{-p) or, equivalently, Pk,ab{p)- 
We have described in Sec. 11.11 that the scahng behavior of these two-point 
functions at an infrared fixed point may be expressed using the occupation 
number exponent k, dynamical scahng exponent z and anomalous dimension 
•q. Using the 0{N) symmetry for = = their scaling forms reacl§ 



^fc,ab(P°, P) 



gL1(/,p) 



1 



fp' 






k I 



s. 



ab 



9 



k^ 



k 



5a 



where the scaling functions / and g^'^ only depend on the scaling ratios 
of the frequency , the absolute value of spatial momentum |p| and 
the scale k. In thermal equilibrium the fluctuation-dissipation relation 
implies the relation ([6]) or ([7]) for k, z and rj. As an application of the 
nonequilibrium renormalization group method, we will determine in the 
following corresponding relations between these exponents at nonthermal 
fixed points using a nonperturbative large- expansion to NLO. 



3.2 Resummed expansion to next-to-leading order 

The search for infrared scaling solutions of the equation for two-point 
functions (j47j) or ( !57l) requires a simultaneous solution for four-point 
functions ( l48i) or ( l59l) . which itself depend on six-point functions etc. In 
order to truncate this infinite hierarchy of equations we employ a resummed 
1/A^ expansion to NLO. A similar calculation for time evolution equations 

^For = the scaling properties are equivalently described by ([5]). 



19 



can be found in Ref. [B]. We write 



n%i^,y) = rr(^,2/)5a. (69) 

and 

+ '^^(z, y) 5^''+'\z - x) d^'^'^y - w) 5aAd 
+ rt*'^*(y, x)5^''+^\y - z) 6^'+'\x - w) 5aad 
+ 5(^^+1) {x - y) {x - z) 5^'^+^) (a; - w) 

X {Sabred + ^achd + (^fccf^ad) • (70) 

The other two- and four-derivatives are written accordingly. For a compact 
diagrammatic notation we will draw propagators along the lines of Sec. 12.31 
and four-vertex functions as 

etc. Plugging (l69l) and (!70|) into the flow equations for the self-energies ( 166!) 
and (|67j) the contributions can be classified according to their scaling with 
A^. The leading contributions ~ turn out to vanish identically. For (l66il 
this is a consequence of the fact that with the vertices fl6Tl) and fl62|) any local 
("tadpole-type") contribution to S-^ vanishes because of /^(^^'^(x, x) = 
and G^{x,y)G^{x,y) = 0. For ([67D the same type of contributions 
cancel because of simple subtraction. Diagrammatically, one is left with 
the equations 




We emphasize that in ( 172|) we wrote the integrated form of the flow equations, 
which "dropped" the scale derivatives appearing in (166!) and (167|) . This trivial 
integration is possible because they represent total derivatives to this order 
in the 1/N expansion as is explained below. 
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For the following it is convenient to introduce one-loop self-energies with 
dressed propagators given bjlll 

(73) 



nr(p) = ^ 



Fk{p - q)Fk{q) - ^pk{p - q)pk{q) 



and a momentum-dependent "effective coupling" \cs,k{p) with 
Acff,fc(p) ^ 1 

A - [i + n^(p)][i + nA(p)] 

= i-n^(p)-n^(p) + (n^(p))V(n^(p))%.... (74) 



For the second line in f[74l) we have used that powers of mixed 
retarded/advanced products I\^{p)Ii^{p) vanish identically. 

The six-point vertex contributes to the flow of F^!^^ only at subleading 
order in the expansion. The leading contribution to the flow described by 
fl59|) with fl70|l generates an infinite chain of "bubble" diagrams with dressed 
propagators when iterated about the short- distance vertices fl6Tl) and fl62|) . 
The fc-integral of this set of flow equations then contains all bubble-chain 
contributions and is trivially obtained from flS^ as 



-^-{>0<->0< + >0<} 

+^Hiv)^{>ao<+>(»;+>ocx 
+ >oo< + >oo< + >oo< } + ••• 

i-nf -n^ + (nf)^ + (n^)^ + ...|nf 



A 

3iV 



Aeff,fc -j-j-iT 



SAT ■ (75) 



^We note that in Fourier space Fu {p) is real while pk [p) is purely imaginary. 
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The second line is derived from the first one in fl75p by simultaneously 
replacing the dressed vertices by their iterated expression to infinite powers 
in Y'^'^'^ and Y'^'^'^ defined in flGlj) and fl62|) . respectively. The first line in the 
above equation depends also on Yff'''^'^{p) and Yff'''^'^{p), as well as Tff'''^'^{p) 
and Vf,'^'^'^ (p) . The first two are given by 

X = X -^^ >OK 

= rr<^{i-n^ + {utr + ...} 

= -3iV + W"^ (76) 

and 

= X -^^ 

= X +(-^'^)';00< +••• 

= rr^{i-n^ + {Y\tf + ...} 

" "l2iV + W^^'^- (77) 

For the last line in (176!) and (177|) we use again that products of retarded 
and advanced self-energies at same momentum vanish. Similar to the short- 
distance vertices (!6T|) and fl62l) . the proper vertex (!76|) plays the role of a 
classical and (177j) that of a quantum vertex. Accordingly, one observes that 
both are related by 

Yf^\p) = i Yt^'^'^ip) , (78) 

and the proper vertex (1771) is absent in a classical-statistical field theory. The 
same holds for the corresponding retarded functions 

rr'\p) = - + nf (p) = I Y^^^ip) . (79) 

With these results the self-energies f l72]) at NLO read 

^kip) = y"Aeff,fc(j)-g)|nf(p-g)F,(g)-in^(p-g)pfc(g)|, 
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n{p) = -^J^\esAp-Q){^kiP-q)Pkiq) + mP-Q)Fkiq)}. (80) 

We observe that this agrees for = with the results obtained from a 
expansion of the 2PI effective action to NLO [22l [H [6] . 

3.3 Fixed points 

We are looking for infrared scaling solutions (168!) of the stationarity equation 
( I65ll . using the approximation (IHOl) for the self-energies. Plugging into (173il 
and (IHOll the scaling form of the propagators ( |68ll . each contribution to the 
self-energies can be classified according to their scaling with powers of k. For 
2 + K > and k > —r] the leading infrared behavior of (p) and (p) is 
described by 

U^ip) ^ ^jF,{p-q)F,{q) (81) 

and 

(p) ^ ^^efr,fc(p - g)nr (p - g)F,(g) . (82) 

The spectral self-energies n^(p) and Sfc(p) are still given by fl73l) and 
flHOj) . respectively, with (p) replaced by flSTj) . We note that the above 
approximation corresponds to the classical-statistical field theory limit of 
the underlying quantum theory, which is obtained by neglecting the quantum 
vertex ( 162|) [251 [T7] . Therefore, we recover the fact that the leading infrared 
behavior of the scaling solutions is the same for the quantum and classical 
theory. Consequently, both the quantum and classical theory belong to a 
universality class characterized by the same critical exponents, which is well- 
established in thermal equilibrium. 

Using (ISTI) and (152]) the stationarity equation (1^ can be written as 

I 5('^+i)(p + l-q- r)X,s,k{p + I) {[Fk{p)pk{l) + Pk{p)Fk{l)] Fk{q)Fk{r) 

J qlr 

-Fu{p)Fk{l) [Fk{q)pk{r) + Pk{q)Fk{r)]} = , (83) 

where we dropped an irrelevant overall factor — A/(18A^). Thermal 
equilibrium trivially solves ([HSD by virtue of the fluctuation dissipation 
relation ([3]). For scaling solutions this leads to the relation ([H]) or ([7]) between 
the exponents z and rj as discussed in Sec. 11.11 It is less straightforward 
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to determine the corresponding relations for nonthermal scaling solutions. 
Following Ref. pQ, this may be efficiently achieved by integrating fl83l) over 
external spatial momentum p and suitable coordinate transformations. In 
this way the problem is reduced to simple algebraic conditions for exponents. 

The integrand of fl83l) contains a sum of four terms, which can be mapped 
onto each other up to proportionality factors using the scaling properties fl^S]) . 
To map the second term in the integrand sum onto the first we employ for 
the frequencies the transformation 

^^4) 

which leads to the absolute value of the Jacobian determinant (p"//*^)^. For 
the spatial momenta we employ 

(85) 

Similarly, for the third term in the integrand sum we take f —* {p^/r^)q^, 
— > {jp/r^)p^, r° {p^/r'^)l^ and p° {p^/r^)r^. Finally, the fourth 
term is transformed with P {p^/q'^)r^, {p^/(f)^^i ~^ {p^/q^)p'^ 

and p^ — >■ {p^/q^) q^- Corresponding transformations are done for the spatial 
momenta analogous to fl85l) . 

It remains to establish the scaling properties of Xes,k{p) defined in fl7^ . 
The propagators fIBS]) determine the behavior of the retarded/advanced self- 
energies defined in (1751) : 

n^'^(/,p) = fc-^.r«'^(|^,M^ , (86) 

where 

A = 4:-d + K-z-r], (87) 
with scaling functions 7i^'^. The infrared behavior of (1741) is then for 

A>0: Aeff,fc(p°,p) = fc^^/(^|^,M^ , (88) 

with scaling function /. For A < the effective coupling becomes trivial with 
AefT,fc(p) — A, on which we comment below. 
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Employing the above transformations to (155]) , and the scaling behavior for 
propagators (1^ and effective vertex (155]) . one obtains for physical correlation 
functions at k = 0: 

X Aefr,fc=o(p + l)Fk=o{p)pk=oil)Fk=oiq)Fk=oir) , (89) 

where 

A = 1 + '- . (90) 

z 

Nonthermal solutions of fl89l) are manifest (1) for A = 0, which corresponds 
to (IE]) given in Sec. 11.11 and (II) for A = —1. The second solution reads 

(II) : K = -ri + 2z + d. (91) 

We finally comment on the case A < in ( !87l) for which the effective 
coupling becomes trivial. One can repeat the above steps with Xes,k{p) — A 
and finds potential nonthermal solutions, which we will denote by an index 
(UV), given by 

8 ^(UV) 

(I): «:(uv) ^ rf_«+^(uv)^!^ 

K J 3 ^ ^ 3 ' 

(II) : «:(uv) = + -^(uv) + H _ (92) 

V ; 3 3 3 ^ ^ 

However, plugging these solutions into (!87|) in order to verify the premise 
dHZD yields (I) (4 - 2r^(uv)) /3 < q and (II) (4 - 2r7(uv) + z^vv)^ /3 < 0. 
For typical values of the anomalous dimension and the dynamical scaling 
exponent this is excluded. 

The type of scaling solutions (l92ll can nevertheless play an important 
role for the physics of (wave) turbulence in classical-statistical systems |28] . 
To understand this, we note that Xes,k{p) — A is physically realized for 
sufficiently large momentum p according to (l74l) . The calculation becomes 
then perturbative and ( l92l) corresponds to the weak-coupling result for large 

if quantum corrections are neglected (i.e. T'^'^'^ = in fl62l) ). The 
perturbative values z^^^^ = 1 and ry(^^) = in (|92|) lead, e.g., in = 3 
to the well-known values k^^^^ = 4/3 and 5/3 [28] . 
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It should be stressed that the UV fixed point solutions (1^ are derived 
neglecting quantum corrections. However, these become important in 
the high momentum regime, where occupation numbers are small. From 
the quantum self-energies flHOl) we observe that the presence of quantum 
fluctuations prevents the UV scaling solutions fl92|) . In contrast, as discussed 
above quantum corrections do not alter the infrared fixed points ([H]) and (1911) 
associated to critical phenomena far from equilibrium. 

4 Conclusions 

Our study reveals a hierarchy of fixed point solutions, which may be classified 
by the exponent k according to (P), ([Z!) and ([H]), or fPT]) . The increasing 
complexity of k from vacuum and thermal equilibrium to nonequilibrium may 
be illustrated in terms of properties of occupation numbers. Vacuum carries 
vanishing occupation number and k is trivially determined by the anomalous 
dimension rj. Thermal equilibrium is characterized by a Bose-Einstein 
distributed occupation number. This distribution is a function of frequency 
and not spatial momenta. As a consequence, in thermal equilibrium k 
depends in addition to t] also on the dynamical scaling exponent z, which 
characterizes scaling with respect to frequency. Nonthermal fixed points 
constitute the least constrained type of scaling solutions in this hierarchy. 
Nonequilibrium occupation numbers depend in general on frequencies as well 
as spatial momenta. We find that k indeed exhibits an explicit dependence on 
the number of spatial dimensions (i, in addition to its possible dependencies 
via 7] and z. 

These results are obtained using the functional renormalization group 
on a closed real-time path. This approach does not build in a fiuctuation- 
dissipation relation, in contrast to Euclidean or complex-time formulations 
at nonzero temperature. As a consequence, it provides a common framework 
to address thermal as well as nonthermal fixed points. Though the derivation 
is in principle a straightforward generalization of standard Euclidean 
formulations, there are important differences which required a thorough 
discussion of the functional aspects of this approach. This includes basic 
relations such as (!33|) . which are crucial for practical calculations of exact 
equations for proper vertices by functional differentiation. 

Nonthermal fixed points represent basic properties of quantum field 
theories, just as thermal or vacuum fixed points, which provide an important 
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means of classification. Further studies are required to reach a similar 
level of understanding as in equilibrium. This concerns, in particular, 
stability properties, which go beyond this work. It is remarkable that 
nonthermal infrared fixed points seem to be approached from a substantial 
class of nonequilibrium instability initial conditions without fine-tuning of 
parameters as reported in Ref. pQ. These results indicate that due to critical 
slowing down the system stays for a long time in the vicinity of a nonthermal 
fixed point, acquiring its properties. Depending on how closely the fixed point 
is approached, subsequent deviations after sufficiently long times finally drive 
the system towards thermal equilibrium. In principle, one may directly tune 
to a nonthermal fixed point without ever deviating from it, which may be of 
measure zero. 

It would also be very interesting to further probe the predicted depen- 
dence of nonthermal fixed points on the dynamical exponent z and on the 
dimensionality of space using lattice simulations in classical-statistical field 
theories in the same universality class. This could be used to verify ([5]) or 
whether more general scaling assumptions have to be made in order to ex- 
plain the observed nonequilibrium properties. Nonrelativistic systems tend 
to have larger z than relativistic theories, which allows a simple variation of 
that parameter as well. Accordingly, a very fruitful direction could be the 
experimental realization of nonthermal fixed point dynamics in cold atomic 
gases, following a parametric resonance or spinodal nonequilibrium instabil- 
ity. 

We thank Thomas Gasenzer, Alexander Rothkopf, and Jonas Schmidt for 
collaboration on related work, and Jan Pawlowski for very helpful discussions. 
This work is supported in part by the BMBF grant 06DA267, and by the 
DFG under contract SFB634. Part of this work was done during participation 
in the programme on Nonequilibrium Dynamics in Particle Physics and 
Cosmology at the Kavli Institute for Theoretical Physics in Santa Barbara, 
supported by the National Science Foundation under grant PHY05-51164. 

A Functional form of propagators 

The identities (l33ll relate second functional derivatives of Wk [J, J\ = 
[J, J, i?^'^''^'^] to those of rfc[0, 0] = r[0, 0; i?^'^'''^'^]. In the presence 
of the i?^'^'^'^ source terms the relations between the fields and J = Jfc[</', (^] 
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and J = Jk[4>,(p] are scale dependent. With (1^ and follows 




6Jk{x) _ r 1 5Jk{x) 6(j){z) ^ 6Jk{x) 54){z) 



) 



which reads in compact matrix notation using fl2^ : 




(93) 



(rr + Rf) wi' + (rf + ) Wi!' = - 



(94) 



Similarly, the other combinations of functional derivatives with J and J can 
be written as 



The above equations can be solved for , i ™^ terms 

of the second derivatives of and R^''^'^'^ , which yields ( !33|l using the 
definitions ([22]). 

B Flow equations for a sharp cutoff 

The flow equations in Sec. 12.31 may be considerably simplified using a sharp 
cutoff function R^{p), which diverges for p"^ < k"^ and vanishes for momenta 
larger than fc^. In Fourier space, the propagators fH2l) can then be written as 



Flow equations for two-point functions given by (H7|) . or fl49p . contain in 
matrix notation 




(95) 



(96) 



(97) 
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Using (pSj) the last expression of (1^ can be written as 

-2fc= A |e(p= _ (rr)"} - e(/ - (r«)" ff* (r 



-1 



= 2kH{p'' - (Tl^j . (100) 
Similarly, 

GtRtGt = 2e5ip'-k')(Tt^y\ (101) 

In the flow equations for two-point functions also appears G^R^iFk = 

G^RfGfvff'G^ and one may be tempted to evaluate this accordingly using 
f llOOj) . However, with (p8|) this naive replacement would lead to ill-defined 
expressions ~ 

5(p2 - P)Q(p2 - k^). Instead, a proper treatment takes into 
account that only the following well-defined sum appears: 

GfRfiFk + iFkR^G^ 



+ i^t + R^y t^i^ + Rty' 



+ (vf + Rf) 'rffrf + Rt) 'tt^Tf + R^' 



In order to integrate flow equations for four-point functions, we employ 

(GMGf) {p - q) vt'Gliq) Vf + G^ip - q) (cf ) (g) ri^^ 
= -2fc2 [5 {{p - qf - k^) e(g2 - k^) + 5{q^ - k^) 6 {{p - qf - k^) } 

X {rry\p-q) rr^(rr)"(.) 

= -2k^ {6 {{p - qf - e( V + 2pg) + 5{q^ - k^) Q{p^ - 2pq)] 

X (rr)~\i>-?)rf (rr)~\?)rf^- (los) 
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Similarly, 

(GfR^Gf) ip - q) T^^^Gtiq) + G^ip - q) vf [GtRtGt) (g) r^^) 

= -2P 1^ _ _ Q(^2 _ ^ ^(^2 _ ^2) Q _ ^)2 _ I 

X (rf)"\p-g)rf (rf)"\g)rl^^ (io4) 

Expressions for four-point flow equations involving iF read 

[gIR^g{) ip - q) T^Miq) rf + Gfip - q) vf (c^i^^F,) (g) vf 

= 2k^5 {{p - g)2 - P) e(g2 - A;^) + _ A;^) 6 ((p - qf - k') } 

X (rr)'\p - q) rl') ((r^)"' (rr)"' ) (g) rl^), (105) 

and 

{GtRtGt) iP - q) iF,{q) + G^ip - q) vf [g^R^iF,) (g) ri^^ 
+ Gt{p-q)Tf (tF.RtGt) {q)Ti'^ 

= 2^ {5 {{p - qf - e(g^ - + 5{q^ - k^) 6 {{p - qf - } 



as well as 

(gM^F, + zF.RtGt) {p-q)Ti'\F,{q)T^^^ 

+ iF,{p - q) (gIrIiF^ + iF^R^Gt) (g) vf 

= -2P 1^ _ ^)2 _ ^2^) Q(^2 _ ^ ^(^2 _ _ ^)2 _ I 

X ((rr)"rr (rr)") (P- g) ((rr)"rf (rr)") (g) rl^) 

(107) 

As exemplified in the last equation of fllOSp . the two theta functions 
contributing to the flow equations for F^'^'* are both independent of k and 
have the same argument but with opposite sign such that they sum up to 
unity after /c-integration. 
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